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Abstract 

We present a canonical quantization of macroscopic electrodynamics. The results apply to inho- 
mogeneous media with a broad class of linear magneto-electric responses which are consistent with 
the Kramers-Kronig and Onsager relations. Through its ability to accommodate strong dispersion 
and loss, our theory provides a rigorous foundation for the study of quantum optical processes in 
structures incorporating metamaterials, provided these may be modeled as magneto-electric me- 
dia. Previous canonical treatments of dielectric and magneto-dielectric media have expressed the 
electromagnetic field operators in either a Green function or mode expansion representation. Here 
we present our results in the mode expansion picture with a view to applications in guided wave 
and cavity quantum optics. 

PACS numbers: 42.50.Nn, 71.36. +c, 81.05.Xj 
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I. INTRODUCTION 



The power of the classical theory of the electrodynamics of continuous media depends 
on capturing the detailed properties of the medium by a small number of spatially-averaged 
fields and effective response functions, such as the linear and nonlinear electric susceptibility. 
The response functions may take a variety of forms depending on the type of system under 
study (e.g., dielectric, magnetic, magneto-electric, optically active, etc.) but in all cases we 
replace the microscopic interactions of an enormous number of charges by a few effective 
macroscopic functions satisfying some general restrictions including the Kramers-Kronig and 
Onsager relations. In this way, many-body systems that would be impossibly difficult to 
analyze directly become easily tractable. 

The use of effective fields and response functions is just as helpful in many types of 
many-body quantum theories and the concept of a quantized theory of macroscopic electro- 
dynamics has held appeal for many authors. However, handling effective response functions 
in quantum mechanics can be challenging because, at least in a unitary evolution picture 
governed by a Hamiltonian, a quantum treatment involving the electromagnetic field is in- 
compatible with dissipation, while dissipation is typically one of the key effects in complex 
systems. In the electrodynamics of macroscopic media, the relation between dispersion and 
dissipation, or loss, through the Kramers-Kronig relations is of central importance, and a 
fully quantum theory must account for it correctly. Moreover, since many important materi- 
als show very strong dispersion, approximate treatments of dispersion can have only limited 
validity. 

Consequently, although quantization of the vacuum field was achieved soon after the 
formulation of quantum mechanics and the corresponding treatment for electromagnetic 
materials was considered soon afterwards, complete formulations of such a theory have only 
emerged in the last decade or so. These theories are timely since metamaterials (MM's) with 
unusual dispersion and significant loss are becoming increasingly common and useful.^ The 
recent studies of spontaneous emission and other phenomena in hyperbolic media [1], for 
instance, suggest that quantum descriptions of MM's and negative media will grow rapidly 
in importance. 

A satisfactory quantization of electrodynamics should have the following properties: it 
should be consistent with the Kramers-Kronig and Onsager restrictions on the response func- 
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tions; recover the Maxwell equations in the classical limit; preserve the correct commutation 
relations between the electromagnetic field operators; and as far as possible accommodate 
otherwise arbitrary constitutive relations. A canonical quantization of the classical theory 
with given constitutive relations meets these requirements. Specifically, one must identify 
a Hamiltonian equal to the energy of the system and expressed in terms of the conjugate 
variables which, when combined with the commutation relations, yields the quantum ana- 
logue of the classical equations; in this case, the Maxwell equations. Efforts to perform such 
a quantization chart a long history commencing with the work of Jauch and Watson on 
a covariant quantum theory of linear, homogeneous, nondispersive dielectrics ^j. Towards 
the same end, Drummond presented a canonical treatment of a nonlinear, dispersive, but 
nonabsorbing dielectric by assuming that the linear susceptibility may be approximated by 



a truncated Taylor expansion over a narrow bandwidth [3]. In attempts to treat causal, 
absorptive media, the main challenge has been reconciling the temporal nonlocality inherent 
to a causal theory of electrodynamics with what must be a temporally local Hamiltonian for- 
malism. This task was completed for a linear, homogeneous, absorbing dielectric by Huttner 
and Barnett 4J who added additional degrees of freedom to the system. In their model, fol- 
lowing the tradition of Hopfield |5|], the electromagnetic field is coupled to a uniform spatial 
distribution of simple harmonic oscillators. In addition, Huttner and Barnett introduced a 
reservoir of oscillators, coupled to the medium oscillators, to facilitate dissipation. However, 
their representation of the canonical variables using spatial Fourier transforms leads to a 
cumbersome theory when applied to inhomogeneous media. Furthermore, the two tiered 



system of oscillators has subsequently been found to be unnecessary jfj]. 

An alternative phenomenological approach to quantization focuses on preserving the com- 
mutation relations of the electromagnetic field operators. This is accomplished by invoking 
the fluctuation-dissipation theorem to introduce source terms into the Maxwell equations 



corresponding to quantum noise currents [7 



t|. In transferring to an operator formalism, 



these noise currents are associated with a set of bosonic fields which ensure that the commu- 
tation relations for the electromagnetic field variables are satisfied. These source terms then 
lead to a Green function representation for the electromagnetic field operators. Several of 
the canonical treatments of macroscopic quantum electrodynamics have been motivated, at 
least in part, by a desire to validate this phenomenological approach. Towards this purpose, 



Suttorp and Wubs 



10j presented a canonical quantization scheme rooted in the Huttner and 



Barnett model, but extended to the case of an inhomogeneous, absorptive dielectric. It was 
then shown by Bhat and Sipe (| that the Huttner and Barnett approach could be refined by 
discarding the medium oscillators and coupling the reservoir directly to the electromagnetic 
field. In addition, magneto- dielectric media have been considered, where the response is de- 
scribed by an electric permittivity e and a magnetic permeability /x. Specifically, a canonical 
treament rjj has made contact with the application of the noise current formalism to a 



magneto-die 
canonically 



ectric 
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121 ]. while spatially dispersive dielectrics likewise have been approached 



Once a quantum theory of macroscopic electrodynamics has been established via the 
canonical route, the door is then opened to the rigorous treatment of quantum electrody- 
namical processes involving dispersive and lossy bulk media. An example of such a process 



is the Casimir-Lifshitz effect 



151 ] whereby forces on solid bodies arise as a result of 



intrinsically quantum mechanical fluctuations in the electromagnetic field. Although the 
best known prediction of this theory is an attractive force between two parallel conducting 
plates jl^j], the possibility of repulsive Casimir forces jl6] has arisen with the consideration 
of left-handed media (LHM). Furthermore, with regard to the study of spontaneous pro- 
cesses, LHM offer novel opportunities in the tailoring of spontaneous emission by atoms 
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18] as well as phase matching in nonlinear optical processes 19l-l2l|. In the absence of 



such materials in nature, LHM are realised through artificial MM's consisting of structures 
engineered, for optical wavelengths, on the nanoscale J22]. In order to treat these materials 
as bulk constituents in an optical system, some process of homogenization must be per- 
formed whereby the electromagnetic response of the sub-wavelength structure is expressed 
by effective parameters for an equivalent continuous medium. Standard parameter retrieval 
techniques, however, frequently return results which appear to violate basic considerations 



such as causality and energy conservation 
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241 ]. initiating a debate on the thermodynamic 



validity and physical interpretation of the effective constitutive parameters [25j . It is argued 



26l | that the problem lies in a neglect of spatial dispersion which may be remedied to some 



extent by modeling the MM as a magneto- electric medium where the polarization and mag- 
netization each depend upon both the electric and magnetic fields. Thus, a magneto-electric 
response may be a general property of any plausible MM realisation of a LHM. Further- 
more, the unusual interaction of light with left-handed MM's arises from their resonant 



properties 



22J, which implies the presence of strong dispersion in the frequency ranges of 



interest for any application which exploits attributes unique to these media. Strong disper- 
sion and, through the Kramers-Kronig relations, strong absorption, is therefore inherent to 
MM based realizations of LHM. A quantum treatment of the electromagnetic field in such 
materials must therefore include causal, magneto-electric constitutive relations representing 
a complete description of dispersion and loss, as opposed to a perturbative approach. 

To our knowledge, a canonical quantization of electrodynamics in a magneto-electric 
medium has not been presented. To achieve this, we must identify a Hamiltonian operator 
which is consistent with macroscopic electrodynamics and yields the desired causal consti- 
tutive relations. As a first step towards such a goal, the oscillator model employed in the 
■nagneto-dielectric case by Philbi, Q w„a,ized in a La grangl an picture to encompass 
a magneto-electric medium by Horsley |27J . However, no construction of the corresponding 
Hamiltonian was attempted. 

In this paper we present a Hamiltonian operator which, with the standard commutation 
relations, is consistent with macroscopic electrodynamics in a causal, linear, inhomogeneous, 
magneto-electric medium. In order to allow for dissipation, degrees of freedom correspond- 
ing to the medium are introduced as bosonic excitations which are then coupled to the 
electromagnetic field variables in a bilinear fashion. This generalization of the interaction 
Hamiltonian constructed previously for a dielectric medium [6] allows for the treatment of 
magneto-electric responses. It is possible to define this coupling to permit a description of 
materials where absorption is absent below a cut-off frequency fl c , such as below the band- 
gap in semi-conductors [6] . In the interests of simplicity we do not make this provision here, 
and rather assume the presence of absorption at all frequencies. Our aim is to obtain the 
dressed eigen-operators of the system, the polariton operators, from which the electromag- 
netic field operators may then be constructed. A significant simplification in the dynamics 
of the system is thus obtained due to the harmonic time dependence of the polariton opera- 
tors. In pursuing this aim we introduce modal polariton operators (i.e. polariton operators 
independent of field point r) and derive corresponding mode field distributions, thus separat- 
ing the spatial dependence of the fields from the time dependent polariton operators. This 
allows us to express the electromagnetic field operators in the form of a modal expansion 
where the polariton operators appear in the place of mode amplitudes. The results thus 
obtained constitute a complete description of the electromagnetic field operators in a broad 
class of a linear, causal magneto-electric media in the absence of a band-gap. 



This paper is structured as follows. In Section III Al we outline the classical theory of 
macroscopic electrodynamics to which the quantum theory must correspond, with particular 
attention given to energy transfer and the general properties of the susceptibility tensors 
describing the response of a causal magneto-electric medium. In Section III Bl we introduce 
our Hamiltonian operator consisting of the electromagnetic field coupled to the bosonic 
excitations of a model medium. We then show in Section III CI how this Hamitonian leads to 
the quantum analogue of the classical theory with constitutive relations expressed in terms 
of susceptibilities which possess the Kramers- Kronig and Onsager properties required of 
their classical counterparts. In Section UTTl we introduce the eigen-operators of the collective 
Hamiltonian and use them to construct solutions for the electromagnetic field operators. 
Finally, a general discussion is expounded in Section IIVI 

II. QUANTIZATION 

A. Classical macroscopic electrodynamics 

1. Independent field variables 

We begin by identifying the key results of the classical field theory which serve as the 
starting point of our canonical quantization scheme. The electrodynamics of continuous 
media in the absence of free charges is governed by the source-free macroscopic Maxwell 
equations, written here in Heaviside-Lorentz (H-L) units, 

D = cVxH, B = -cVxE, (la) 
V-D = 0, V-B = 0. (lb) 

The field variables D(r, t), B(r, t), E(r,t) and H(r,t) are the electric induction, magnetic 
induction, electric field and magnetic field, respectively. A dot above a quantity denotes a 
time derivative, and c is the speed of light in vacuo. Conversion to SI units is effected by 
replacement of D, B, E, and H, by D/y/e~o, B/y^/zo, y/eo E, and y/JIo H, respectively, where 
Eq is the permittivity of free space and /io is the permeability of free space. Although the 
choice of H-L units is somewhat unorthodox, it is particularly convenient in the consideration 
of magneto-electric media in that the electromagnetic field variables all share the same 
dimensions. From the outset we may identify fllbl) as initial conditions for D and B since, 
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on account of (flal) . if they are satisfied at one time they are satisfied at all times. In 
addition to (JT|), a complete description of the field dynamics requires that one pair of the 
field variables be treated as independent and the remaining pair be expressed as functions 
of them through a set of constitutive relations. In doing so, there is a freedom regarding 
which two quantities are chosen as the independent pair. Here we take D and B to be the 
independent fields. In order to justify such a choice we write down the standard expression 
for the incremental change in the energy density as a result of changes in the fields alone, 

m = E • dD + H • dB. (2) 

In the absence of dispersion the relationship between the field variables is local in time and 
E and H may then be written as 

which then allows ([T]) to be rewritten solely in terms of D, B, and the energy density 
viz., 

D = CVX 5B' B = " CVX 5D' (4a) 
V ■ D = 0, V ■ B = 0. (4b) 

The pairs (E, B) and (D, H) are commonly viewed as primary and subsidiary variables, 
respectively (e.g., |28|). However, ()2])-(j4|) suggest that, for a Hamiltonian picture of dis- 
persionless macroscopic electrodynamics in the absence of free charges, the natural choice 
of independent field variables is the pair (D,B), as first noted by Born and Infeld j^]. 
Motivated by this result, as well as the advantage of working with transverse fields, we 
extend this choice to the present treatment where dispersion and loss are included, and 
no longer follows directly from fl2]). Nonetheless, we note that the quantum analogue of ([3]) 
holds with the energy density ^ replaced with the Hamiltonian density corresponding to 
the Hamiltonian operator H to be presented in Section III B[ 

The required constitutive relations must therefore express the fields E and H in terms 
of D and B at all times. With the standard definitions of the polarization P(r,t) and 
magnetization M(r, t) as 

P = D-E, M = B — H, (5) 



TABLE I. Definitions of some common classes of media through the dependencies of the polarization 
P and magnetization M upon the electric induction D and magnetic induction B, when the latter 



pair are chosen as the independent fields. 

Medium class P M 

Dielectric P{D} 

Magnetic M{B} 

Magneto-dielectric P{D} M{B} 

Magneto-electric P{D,B} M{D,B} 



this is equivalent to writing 

P = P{D,B}, M = M {D, B} , (6) 



where the temporal non-locality is implied (SI units are obtained by the replacements P — > 
P/a/^o and M — > y^uo M). In the absence of explicit spatial dispersion, the dependencies of 
P and M upon the independent fields D and B in represent the most general case of an 
anisotropic, gyrotropic, magneto-electric medium. The definitions of several other common 
classes of media are summarized in Table HI 

2. Susceptibility tensors 

The choice of independent fields, and therefore the form of the constitutive relations in 
(|6]), have ramifications for how the transfer of energy between the electromagnetic field and 
the medium is viewed. By considering this energy transfer we now derive some properties 
of the medium response functions. Standard manipulation of (TTa]) yields 

- V S EH = E D + H B, (7) 

where Seh = cE x H is the Heaviside-Lorentz form of the electromagnetic energy flux 
associated with E and H. With the present choice of independent variables it is natural to 
define the energy Udb associated with the fields D and B in the volume V, 

U DB = / dV \ (D • D + B • B) . (8) 
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We then insert fl3j) into (J7|) and integrate over all time and over the volume V to obtain 

- f dt <f dA n • S EH = AU DB - / dt f dV (p • D + M • b) , (9) 

J-oo Js J-oo Jv ^ ' 

where AUdb is the total change in the independent field energy, and n is the outward 
unit vector normal to the surface S. There has been much discussion regarding the correct 
definition of the Poynting vector as the representation of electromagnetic flux (e.g., 30] ). 
By taking S to lie exclusively in the vacuum we avoid this controversy in that, for the 
purposes of the integral on the LHS of (0), all the possible flux vectors are equivalent; i.e., 
Seh = Seb = Sdh = Sdb = cD x B. Assuming that the only flux of energy across S is of 
an electromagnetic nature, we may then use fl9]) to write the total change in the energy of 
the system enclosed within V, AU to t, as 

AU to t = - / dt I dA n • S DB = AUdb + / dt f dV (p • D + M • b) . (10) 

J-oo JS J-oo Jv ^ ' 

The second term on the RHS of (fTUj) represents the contribution to AUtot due to the presence 
of a medium within V. In order to relate this change in energy to the response of a magneto- 
electric medium, we write explicit expressions for the constitutive relations in OH]) by defining 
the real susceptibility tensors T au (r,t), with a, v = e,m. These relate P and M at time t 
to D and B at all times through 

/oo 
dt [r cc (t - 1') ■ D(t') + r cm (t - f) ■ B(t')] , (n) 
-oo 

/oo 
dt [r me (t - • D(t') + r mm (t - 1') ■ B(t')] , (12) 
-oo 

with the Fourier domain representation 

P(w) = r co (u;) ■ D(w) + r cm (u;) ■ B(w), (13) 
M(u) = T mc (cu) ■ D(w) + r mm (a;) ■ B(w). (14) 

In order to identify the dissipative and non-dissipative parts of the susceptibility tensors we 
define 

rs P («) = i[i w ( w ) + f w (w)], (is) 
r^) = l[r'>)-f^*H], (16) 

where an over-bar denotes a tensor with Cartesian components obtained from those of 
the unbarred quantity by a matrix transpose (recall that a, v = e, m). Thus an over-bar 



combined with a star indicates the Hermitian transpose (the superscript f is reserved for 
the adjoint of an operator). This allows T au to be written as 



diss 



(17) 



We note that for a = u, (u;) and r^ s (u;) are the Hermitian and anti-Hermitian parts, 
respectively, of T aa {uj). We may then employ f lT3|) - f|T6|) to express the second term on the 
RHS in ([TO]) as 



v 



d* / dV P-D + M-B 



-lUJ) 



dV 



D*(w) B*(w) 



diss 



diss 



D(w) 
BM 



where the Fourier transform of an arbitrary function f(t) is defined as 



dt e luJt f(t). 



, (18) 



(19) 



Now consider a transient interaction of the electromagnetic field with a medium of finite 
extent, such that S may be assumed to lie in the vacuum and AUdb — (Udb(^) — for 
t = ±oo, say). From (fit)]) we see that the RHS of ([IB]) is unambiguously the total energy 
transferred between the electromagnetic field and the medium. Thus, only the medium 
response described by r^ s (o;) leads to gain or dissipation of electromagnetic energy, and 
T^p(u) describes the non-dissipative interaction. 

Further symmetries of the susceptibility tensors may be determined through applying 
time reversal. We impose the following properties upon the fields in the time domain: 



B, 



{D}_.^ = D, {B}_^ 
{P}_^ = P, {M}_^ = -M, (20) 

where {}_^ represents the operation of time-reversal upon the enclosed expression. Addi- 
tionally, the application of time-reversal in the Fourier domain corresponds to the replace- 
ment u — > —uj, as well as the appropriate transformation of any parameters; e.g., in the case 
where the susceptibility tensors are dependent upon an ambient magnetic field B we have 



{IX(w,Bo)}_^ = {rSL(-w,B )}_ 



dissv 



-uj, 



-Bq). Using {}_b to represent the 



time reversal of all such parameters, we may then use f fl8|) . along with the reality condition 
for the fields, to obtain 



{r3L(w)}_ 



diss 



{iX»Lbc 
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diss 



a v. 



(21) 



Expressions analogous to (121]) for TZ sp (ou) are obtained by imposing causality, 
to setting T au (i) = for t < 0, which results in the Kramers-Kronig relations 



'his amounts 



31 



32| 



-oo 2vr w- 



rZM = 2 ^l ^-^r> ( 22 ) 



TZM = ^ ^z^—r^ (23) 

where & indicates a Cauchy principal value integral. Combining ( )2T|) - (l23|) then leads to 

{r- p (o;)}_ Bo = f Sp(w), {nr sp (o;)}_ Bo = -f Sp( W ), a + v. (24) 
From (II 7p . (I2ip . and ( 124)) we may now identify the Onsager relations 

{r^l-Bo = f ^ {TU = -f v \ (25) 

which hold in both the time and frequency domains. It should be noted that in deriv- 
ing the Kramers-Kronig and Onsager relations expressed in f l2"2"j) . f l2"3"j) . and (l25p. only the 
macroscopic Maxwell equations, the behaviour of the fields under time reversal, and the 
assumption of causality have been used. Therefore, these conditions upon the susceptibility 
tensors represent fundamental properties of a causal, electromagnetic medium, and they 
must be reflected in a valid quantum theory. 



B. Hamiltonian operator 

We now turn to the construction of a Hamiltonian operator that, with the standard com- 
mutation relations, leads to the quantum analogue of (CEJ) with constitutive relations of the 
form given in ( ITTj) and ( IT2l) . The associated susceptibility tensors must satisfy standard 



Kramers-Kronig and Onsager relations. Some Hami 



electrodynamics have proceeded from a Lagrangian (4, |10|, |ll|. However, it is sufficient to 
provide a Hamiltonian operator directly. Indeed, transforming from the Lagrangian to a 
Hamiltonian leads to a complicated field-medium interaction that is difficult to diagonal- 
ize. Instead we directly construct a Hamiltonian which describes a rather general class of 
magneto-electric media. This allows us to define the various couplings straightforwardly in 
terms of the canonical variables, and avoid the complications of transitioning from a La- 
grangian to a Hamiltonian picture. In addition, the classical value of the Hamiltonian must 
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tonian formulations of macroscopic 



be equal to the energy of the system, which is satisfied in our theory below by construc- 
tion. The full system consists of two linearly coupled subsystems representing the vacuum 
electromagnetic field and a medium. The corresponding Hamiltonian is therefore of the form 



H = H emf + H med + H 



mt ? 



(26) 



where H em f, H me d, and Hj nt are the electromagnetic field, medium, and interaction Hamil- 
tonians, respectively, and a hat denotes an operator. In what follows, all operators may 
be presumed to commute unless specified otherwise. In the Heisenberg picture, the time 
evolution of an arbitrary operator O is governed by the equation 



ihO 



0,H 



(27) 



where [ , ] represents a commutator, and the dot denotes a total derivative with respect to 
time (in that the components of the field point r are not dynamical variables). 
The Hamiltonian for the electromagnetic field is 29j (c.f. (JH])) 



H 



cmf 



dV 



D(r,t) -D(r,t) + B(r,t) -B(r,t) 



(28) 



and the equal time commutation relations (ETCRs) for the components of the field operators 



are 
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d 

A(r, t), Bj(r', t) = ihce ikj —5(r - r ) 



(29) 



where Cartesian vector and tensor components are indexed by Latin subscripts, a sum over 
repeated indices is implied, e is the Levi-Civita pseudotensor, and 5(r) is the Dirac delta 
function. 

The construction of the medium and interaction Hamiltonians requires some discussion. 
Previous models for a medium which exhibits both an electric and magnetic response have 
typically involved independent electric and magnetic subsystems. In the application of the 
phenomenological approach to a magneto-dielectric medium [jjl , the noise polarization and 
magnetization are made to originate from independent bosonic vector and pseudovector 
fields, respectively. Such a separation of the medium is deemed appropriate in modeling 
materials where the electric and magnetic responses arise from physically distinct material 
constituents or degrees of freedom. Similarly, in the canonical treatment of a magneto- 
dielectric 11] this prescription is reflected in the introduction of two separate sets of harmonic 



oscillator fields in the model for the medium. Following on from this work, the identification 
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of these fields as corresponding to electric and magnetic oscillators was then made explicitly 
in the construction of a Lagrangian for a magneto-electric medium [23]. This was effected 
by associating the symmetry properties of the oscillator amplitudes under spatial inversion 
and time reversal with those of the electric or magnetic field, as appropriate. We have found 
that a magneto-electric response may be obtained with a single set of oscillators. However, 
such a response is severely restricted (e.g., the cross-coupling is completely fixed by the 
dielectric and magnetic responses in the isotropic case). Here we present the most general 
form of magneto-electric response obtainable within the context of the established harmonic 
oscillator model. Since the coupling of the field to multiple continua is already implied in the 
vector nature of the oscillators, it is a straightforward matter to add additional degrees of 
freedom to the medium. We pursue this by including an arbitrary number of such oscillator 
sets (to represent a MM with electromagnetic resonances associated with several material 
constituents, for instance). In general, we may expect each class of oscillators to exhibit 
a magneto-electric coupling, with no purely electric or magnetic character, and we model 
the medium with N sets of vector operators h\Q(r,t), A = 1,...,N, representing bosonic 
excitations associated with each frequency fl > and field point r. The required symmetry 
properties of the response are then imposed on the coupling coefficients between the field 
and medium. The corresponding Hamiltonian for the medium is 

/r°° HQ 
dV — mb{ Q (r,t)-b X a(r,t), (30) 
A JO 71 

and the vector components of the medium operators obey the ETCR 

bx H (r,t),b{ l{Vj (r f ,t)\ = 2ir6xx'S ij 5(Q-Q') 6(r-r'). (31) 

where 5w and 5y are Kronecker deltas. The interaction Hamiltonian is then constructed 
as a spatially and temporally local, bilinear coupling between the vacuum electromagnetic 
field and medium operators, viz., 

/POO JQ ,- -, 
dVD(r,t).^/ — [AUT,n)-b xn (T,t) + A c x *(v,n)-b[ n (v,t) 

/roo HO r 1 
dVB^-^j ^ [^(r,0)-b AC (r,t) + Ar(r,0).b^(r,t)_ , 

(32) 

where the complex valued, second rank proper tensors A^(r, fl) and pseudotensors A™(r, fl) 
are defined for positive fl. Finally, we assume that the microscopic dynamics that underlie 
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our macroscopic picture lead to a Hamiltonian which is symmetric under time reversal. We 
therefore require that A e x — y A e x * and A" 1 — y —A™* under this operation. This requirement 
also ensures consistency with the classical theory (see (157)) and (155)) below). For the alter- 
nate representation where the medium fields are represented as harmonic oscillators with 
coordinate operators q xn (r,t) and conjugate momenta p xn (r,t) (e.g., j(| the form of 
(132)) corresponds to the coupling of D and B to both q xn and p xn . Ultimately, the coupling 
tensors A x and A™ in (1321) are to be determined by the measured or calculated susceptibility 
of the medium. However, making the replacements A x —y A x ■ XJ* X and h\n —y XJ\ ■ h\n, 
where a = e, m and the tensor XJ X represents an arbitrary unitary transformation, leaves the 
Hamiltonian unchanged. This represents an inherent freedom in the model for the medium. 
Thus, inserting (1281) . (I3UI) . and (152")) into (1261) we may write the full Hamiltonian explicitly 

as 



H = - J dV [D(r, t) ■ D(r, t) + B(r, t) ■ B(r, t) 

/r°° c\n 
dV —hnb{ n (r,t).b xn (r,t) 
. */ 



h? / dVD(r,t) • 



dQ 
dtt 



A c A (r,fi)-b An (r,t) + Ar(r,fi)-bi n (r,t) 



hk I dVB{r,t)-J2 /_ ^ A-(r,fi)-b An (r,t) + Ar(r,fi)-bL(r,t) .(33) 



C. Dynamical equations and constitutive relations 



To demonstrate the consistency of the Hamiltonian system presented in Section Til Bl with 
macroscopic electromagnetism, we first insert D and B into ( 1271) and use ( )29l) to obtain 



D = cVxB-cVxfi^/ ^ A^(0).b A n + Ar(^)-bi n 



271 



(34) 



and 



B = -cV x D + cV x hh Y 



dn 



A e x (n)-b xn + A c x *(n)-b 



xn 



Recalling (J3J) and making the operator definitions analogous to 



(35) 



P = D-E, M = B — H, 



(36) 
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we are led to identify the polarization P(r, t) and magnetization M(r, t) of the medium as 



A 







2tt 



A e A (r, O) • b AQ (r, t) + A C A * (r, Q) • b AQ (r, t) 



(37) 
(38) 



Note that consistency of ( 157)) and (1551) with (120)) follows from the time reversal properties 
imposed upon the coupling tensors. Inserting (156)) - (138)1 into (151)1 and (135)) we obtain the 
quantum analogue of the Maxwell curl equations, ( ITa)) . as 



D = c V x H, 



B 



-c V x E. 



Additionally, the conditions 



V-D = 0, 



V-B = 0, 



(39) 



(40) 



must be enforced independently at this stage to reproduce the full set of macroscopic Maxwell 
equations. The description of the dynamics of the electromagnetic field operators contained 
within ( 156l) -(l40l) remains incomplete, however, until we establish constitutive relations of 
the form (ITU and (TT2|) . 

To identify the required constitutive relations for the medium we must determine the 
contributions to the dynamics of b A f}(r,t) driven by the electromagnetic field. Having thus 
expressed the medium operators in terms of D and B, we may substitute the results into 
(157)) and (1581 to yield the field-induced polarization and magnetization. The dynamical 
equation for the medium operators follows from (127)) . (15T1) . and (155)) . as 



b An (r,t) = -iQb xn (r,t) + itr* [A c /(r,ft) . D(r,t) + A^fr fi) ■ B(r,t) 
Integrating this equation directly for some initial time r < t we obtain 

/CO 
dt > Q(t - t') e - m (*-*') [A c A *(r, Q) ■ D(r, if) + A^Cr, fi) • B(r, t') 

+b xn (r,r) e-' 1 ^. 



(41) 



(42) 



where the causality in the relationship between the electromagnetic field operators D and 
B, and the medium operators b X n is clear. The use of the Heaviside step function 6{t) has 
allowed extension of the upper limit of the integral to +oo. 
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Returning to the task of identifying the constitutive relations, we substitute ( 14*2]) in ( 137]) 
and fl38|) to obtain expressions for the polarization and magnetization operators, 



P(r,t) 



M(r,i) 



dt' 



r ee (r, t-t')- D(r, t') + r em (r, t-t')- B(r, t') 



+P (n) (r,t), 



dt' 



r me (r, t-t')- D(r, t') + r mm (r, t-t')- B(r, £') 



+M( n )(r,t), 



(43) 



(44) 



where the r CT ^'s are identified below, and we have defined the noise operators 

"°° dO 



P (n) (r,t) =t£ J2 



x 



o 



2tt 



A A (r,0)-b An (r,r) e"* ^ . 

+A c A *(r, O) • b^r, r) c «<*-> 
A^(r,n).bAn(r,r) e" 1 ^ 

+Ar(r,n).bJ n (r,r) 



(45) 



(46) 



The domains of the integrals in (1431) and (14*41) have been extended to — oo by choosing r 
such that t — r > tr, where tr is the finite response time of the medium. The expressions 
( l4"3"j) and ( |44"|) form the quantum analogue of the constitutive relations (TXTT) and (fl"2l) . as 
desired, with the addition of the noise terms p( n )(r, t) and M^(r,/f:) involving the initial 
conditions for the b A ^(r,t) operators. Such noise operators are a hallmark of dissipative 
quantum systems wherein they act to preserve the commutation relations by compensating 
for the otherwise dissipative decay of the coupled operators. 

The causal susceptibility tensors r CTy (r, t) in fj43l and (j4"4"l) are related to the coupling 
tensors A A (0, r) (a, v = e, m) through 

f°° HO 

r"(r,t)=0(t)V / — [2^{A^(r,0)-Ar(r,0)} sin(Oi) 

- 2^{A A "(r,0) • A A *(r,0)} cos(Ot)] , (47) 

where the symbols &{} and ^{} represent the real and imaginary parts of the enclosed 
expressions, respectively. The Heaviside step function on the RHS of (J4Tjh along with the 
time reversal properties imposed upon the A A 's, ensures that the identification made in 
equating the LHS and RHS of (14T1) is consistent with the requirement that the T au, s satisfy 
the Kramers-Kronig and Onsager relations as expressed in f[22l . f[23|) . and (125]) . In the 
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Fourier domain, r au (r, u) may be separated as in ( 1171) . and by taking the Fourier transform 
of (jUJ) we may identify ^^(r,^) and r^ s (r, u), defined as in ffl5l) and ffl6l) . respectively, 
in terms of the coupling tensors as 



Thus, ( H3|) -( l49|) demonstrate how the temporally local coupling in (132|) is related to the 
causal response of the medium. 

So far we have demonstrated that the Hamiltonian and accompanying ETCR's intro- 
duced in Section III Bl lead to the quantum analogue of the macroscopic Maxwell equations 
with magneto-electric constitutive relations involving susceptibility tensors which obey the 
Kramers-Kronig and Onsager relations. In addition, expressions for the non-classical noise 
operators PW and MW have been obtained, which fully determine their commutation rela- 
tions and time evolution properties through dependence upon the initial conditions for the 
medium operators in the form hxn( r , T ) e~ 4f ^*~ T * ) . Thus, Section HT1 constitutes a canonical 
quantization of macroscopic electrodynamics in a causal magneto-electric medium, within 
the restrictions placed upon the susceptibility tensors by their relationship to the coupling 
tensors in fj47|) -f l49|) (see discussion in Section HV|) . 

III. POLARITON OPERATORS 

The quantization of Section [TT] has provided us with the dynamical equations for the 
electromagnetic field operators, and we are thus in a position to identify their solutions. 
The resulting expressions for D, B, E, and H constitute the main results of the present 
work and are derived below in the form of the mode expansions floTl) . fl62|) . fl86|) and fl87j) . 
with the mode fields following from solutions to ( 180|) . In obtaining these solutions we employ 
standard methods of modern quantum optics, which we now proceed to describe. 

In the Heisenberg picture the dynamics of coupled quantum systems are often simplified 
by determining the eigen-operators of the full Hamiltonian which exhibit harmonic time 
dependence. When associated with dissipative systems, the application of this procedure is 




(48) 



and 




(49) 
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known as Fano theory [34], |35[ . In the present context of macroscopic electrodynamics we 
term the eigen-operators of H as the polariton operators which are required to satisfy the 
eigen-equation 

mo n {t)= [6n(*),H] , (50) 

whence, with (l2"Tj) . On(t) = O^(0) exp(— iQt) follows. However, ([50]) tells us nothing about 
how the spatial degrees of freedom in the Hilbert space implied by (I3"3"]) are to be accounted 
for in the labeling of the polariton operators. 

In the canonical quantization of the vacuum electromagnetic field the concept of photons 
arises in connection with the plane wave modes of the vacuum. These photon modes are 
associated with the purely time dependent eigen-operators of the vacuum field Hamilto- 
nian H em f, which evolve as a n k(^) = a„,k expf— iu^t). This allows the electromagnetic field 
operators to be written as expansions of the form |36( 

D(r, t) = J2 [a«k(f) D uk (r) + a[ k (t) D^(r)] , (51) 

B(r, t) = J2 [a«k(*) B uk (r) + a\^t) B^(r)] , (52) 

where D u k(r) and B u k(r) are the vacuum wave modes with wavevector k and polarization 
index u = 1,2. The discrete sum over k follows from normalization to a finite box with 
periodic boundary conditions. In the extension of this approach to a linear, inhomogeneous, 
but non-dispersive medium (e.g., js]), the forms of ( 15 ip and (152]) are preserved with the plane 
wave solutions being replaced with the electromagnetic modes of the structured medium. 
Such a description is particularly appropriate in the context of guided-wave optics and 
photonics, where the spatial modes of waveguiding structures and optical cavities form the 
natural language for the dynamics of the system. For this purpose it is desirable to maintain 
the modal approach when extending the quantum treatment of the electromagnetic field to 
Kramers-Kronig media, as is done here. We therefore introduce modal polariton operators 
in order to build a description of the electromagnetic field operators D and B analogous to 
(I5T]) and (I52p . The polariton operators are labeled by the indices n and fl, the nature of 
which follows from consideration of the Hilbert space implied by the form of ( ]3"3"]) : the index 
Q is continuous, while the precise nature of the index n, which represents the spatial degrees 
of freedom, is determined by the geometry of the system along with the boundary conditions 
imposed: ultimately we shall find that n labels spatial electromagnetic field distributions 
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corresponding to the polariton modes in analogy with the labels (u, k) for the photon modes. 
For convenience, as in the vacuum mode case, we assume normalization within a finite box 
with periodic boundary conditions and therefore treat n as a discrete index; the general- 
ization of this prescription to multiple or continuous indices is straightforward and may be 
made later as required. 



A. Transverse and longitudinal response polaritons 



From the outset we partition the modal polaritons into two classes: transverse response 
polaritons (TP's) and longitudinal response polaritons (LP's), leading to the Hamiltonian 
form 

fi = E/ or ^ 4,(0 &»(*)+££/ ^^4„W^W, (53) 



n JU X n " 



where cn n (t) and sxn n (t) correspond to the TP and LP operators, respectively. The TP 
operators represent the collective field- medium excitations and form a single class. The LP's 
represent excitations of the medium which do not couple to the transverse electromagnetic 
field and therefore constitute at most N subclasses for each subsystem of the medium. By 
assumption, the polariton operators satisfy the equations 



Mlcn n (t)= Cnn(t),H 
Mlsxsm(t)= s X nn(t),R 

which imply CQ n (t) = c<^ n e~ lClt and S\n n (t) = s X nn e~ tQt . We also impose the ETCR's 



2n5 nn i5{yt — Q ), 

2ll5xXi5nn'5(VL - Q'), 

0. 



(54) 
(55) 

(56) 
(57) 
(58) 



The last of these ETCR's establishes the formal separation of the two classes of polariton 
operators, the classification of which follows from the definition of the LP modes as those 
associated with configurations of the medium which do not interact with the transverse D 
and B fields through the interaction Hamiltonian. This may be expressed in the form of the 
condition 

^ n ,H^l =0, (59) 
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where 



H 



(A) 
hit 



-ft* / dV t>{r,t) 



dQ 
2^ 



A c x (r,Q) ■ b xn (r,t) + AT(r,Q) ■ b{ Q (r,t) 



dVB(r,f). J ^ [AA(r,fi)-b Af ,(r,t) + Ar(r,fi)-bL(r,t)J. (60) 

The partitioning of the polaritons into TP's and LP's combined with the transversality of 
the electromagnetic field operators now allows us to derive a number of useful results. 

Since the polariton operators are the eigen-operators of the full Hamiltonian we may 
write D, B, and b as expansions of the form 

roc dQ 



£>M) = E 



2tt 
dQ 



cn n (t) D nn (r) + cL(t) D^ n (r) 



2tt 

dQ' 
2^ 



B(r,t) = J2 / oT &»(*) B nn (r) + 4,(t) B^(r) 



(61) 
(62) 



cn'n(*) «An(r, Q', Q) + 4 n (t) /3* n (r, tf, ft) 



(63) 



where Do n (r), Bjj n (r), CKAn(r, f2', Q), (3\ n (j, Q', ^), and p\ n (r, Q) are vector field coefficients. 

To obtain the inverse transformations corresponding to (l6"TT)-(l63l) we construct nominal 
expansions of the polariton operators in terms of the subsystem operators D, B, and b. 
Using these expansions to evaluate the ETCR's between the subsystem and polariton op- 
erators (e.g., [D(r), CQ n }) and comparing the results with the same ETCR's evaluated using 
(16"I|)-(163]), we obtain 

= ^ f dV f E nn*(r) • D(r,t) + H^*(r) • B(r,t) 



hQ 



+ Y1 dv 



dQ' 
2^ 



a An (r, Q, Q>) ■ h XQ , (r, t) - /3 An (r, Q, fi') • b{ Q , (r, t) 



W*) = / dVp An (r,fi)-b Af ,(r,t), 



(64) 
(65) 

where, due to the ETCR (I29I) . the vector coefficients E^(r) and H^(r) are related to 
D nn (r) and B nn (r) through 

- iQB Qn (r) = c V x HTP(r), -zQB nn (r) = -c V x E^(r). (66) 

We will eventually solve for the coefficients Djj„(r) and Bn n (r). Thus (1561) only defines the 
transverse parts of E^(r) and H^(r), with the longitudinal parts being unconstrained. 
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Convenient choices for E^(r) and H^(r) are made in (!78|) and ( 179]) below so that they 
correspond to the coefficients of the TP operators in the polariton expansions of E and H. 

Using the expansions ( 16T]) -f l65]) . the conditions fl58|) and f )59|) defining the partitioning 
of the TP's and LP's may be expressed as requirements upon the expansion coefficients. 
Substituting fl6"Tj) - (l63"]) in fl32|) we may reexpress the definition f )59|) of the LP's as 

J dV p* Xn ,(r, ST) ■ [A c A *(r, fi') • D Cn (r) + M) ■ B nn (v)} = 0, (67) 



for all A, n,n', and Q,Q'. Likewise, by substituting ( I64p and (I65p into the condition 
and evaluating the commutator we may reexpress the formal separation of the TP's and the 
LP's as 

J dV p* Xn ,(r, Q') ■ a Xn (r, Q, Of) = 0, (68) 
for all A, n, n' and Q, fl'. Use of ( 165)) in ([37)1 implies 

J dV pl n (T,n).p Xn ,(r,n) = 5 nn , (69) 

Inserting the expansions f )6T]) -f l63|) into the dynamical equation for the medium operators 
( |4T|) and evaluating the ETCR with c^ n and c^ n we obtain 

(n' - n) cx Xn (r, n, n') = ri [A c A *(r, n') ■ D nn (r) + A%*(t, n') ■ B n „(r)] , (70) 
(Q' + Q) (3 Xn (r, Q, Q') = ft" a [A c A (r, Q') ■ D nn (r) + A?(r, Q') ■ B nn {r)] . (71) 



Following the approach of Fano 
fl68|) to obtain the solutions to (170)) and flTTj) as 

, 1 



34j as extended by Bhat and Sipe we may use (I6"7j) and 



0!An(l\ ^ ^'1 = ^ 2 



+ Z Qn 5(0' - Q) 



x [A c A *(r, fi') ■ D nn (r) + Af(r, fl') • B nn (r)] , (72) 
/3 AB (r j n,n / ) = fH— ^ K(r,^) • D nn (r) + A x (v,tt) • B nn (r)] . (73) 

The complex scalar Zn n parametrizes the resonant interaction between the electromagnetic 
field and the medium as a whole and is therefore independent of A. It represents a generalized 



contribution from the pole at Q = Q' originally introduced by Dirac [37J • We find that Z, 



which is undetermined at this point, emerges as the eigenvalue associated with solutions of 
a generalized Hermitian eigenvalue problem. No pole is considered in obtaining (173)) since 
Q and Q' are positive. 
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Thus, the expressions (1721) and (1751) have reduced the unknown quantities associated with 
the expansions (16T1) - (163|) to just the vector fields p An (r, f2), D^ n (r), and Bn n (r), and the 
scalar Zq u . The procedure for determining p An (r, O) was given by Bhat and Sipe for a 
dielectric medium 6|. Although we consider a more general coupling here, similar principles 
apply and solutions for p An (r, fl) may be obtained from consideration of (KT7j) and (1551) . 
along with the transversality of Dn n (r) and B(j n (r) (which is derived below). However, we 
do not consider this solution procedure explicitly here as only the transverse electromagnetic 
field operators correspond to optically measurable variables, and any coupling between the 
electromagnetic field and an atomic system, for instance, is mediated by the same. We 
therefore restrict our attention to the TP operators and the explicit construction of D(j n (r) 
and Bn n (r), as these alone contribute to D and B through (KH| and (1521) . 



B. TP mode fields 



The vector fields D(j„(r) and Bo„(r) are of particular importance as they play the role 
of TP mode fields. From (l6T|) and ( l62l) it may be observed that, since the time evolution 
and commutation properties of the TP operators are known, determination of the TP mode 
fields constitutes a complete description of the electromagnetic field operators D and B. 

To determine the TP mode fields we insert (l6Tl) - (l63l) into (1341) and (155]) and evaluate the 
ETCR of both sides with c^ n to yield 



iQ D n „(r) = c V x B nn (r) 



an' 

o 2tt 



+Ar(r,n').(3 Xn (v,n,Q')} 7 (74) 



iQ B n Jr) = -c V x D n „(r 



+c V x m V / — [A e A (r, Q') ■ a An (r, Q, Of) 



(75) 



Substituting Eqs. ( J72|) and ( J73|) into Eqs. ( J74|) and (175|) . and recalling the definitions in (fTTj) . 
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HJ), and (JUJ), then leads to the TP mode field equations 



- in Dnn(r) = c V x { [l 



rdisp( r > ^ 



B nn (r)-f^;(r,fi)-D nn (r)} 



c^V x [lX(r, O) ■ B^(r) - f e ™(r, fl) . D nn (r 



Z7T 



-in B nn (r) = -c v x { [i - r- p ( r , n)] • n nn - ix P (r, n) • B n „} 

+c4^V x [I- s (r, n) • Dnn(r) + r^(r, fi) • B n „(r)] . 

Z7T 



(76) 
(77) 



If we set = in these equations correspond to the classical, source free Maxwell curl 
equations in the frequency domain, with the associated complex frequency solutions. How- 
ever, this is inconsistent with unitary evolution. Instead we seek complex solutions for Zn n 
which determine an augmented medium response, thus ensuring the reality of the polariton 
frequency Q. From a related perspective, it is shown below in Section IHI CI that the devi- 
ation of Zn n from a value of in determines the contribution of the corresponding polariton 
mode to the noise polarization and magnetization. 

A consequence of (ITS"]) and (ITT)) are the conditions V • D^n(r) = V ■ B A n(r) = 0, from 
which fHOj) follows. Thus, construction of the field operators D and B according to ( IBTj) 
and (I6"2"|) ensures their transversality, and the condition (130)) need no longer be enforced 
explicitily as an initial condition. Comparing ( 166)) with (1T61) and (1771) we are free to identify 



E^(r) 



[1 - IX p (r, Q)] • D An (r) - T^ p (r, fl) • B A , n (r) 
/Ml T^s(r, O) ■ D An (r) + r^(r, O) • B An (r)] , 



— c 



?7T 



H 



TP/ 



[1 " C P (r, «)] • B A n(r) - n-(r, O) ■ D Af ,(r) 
A<, Tr (rj n) _ B ^ (r) + n) . DAQ(r)] _ 



— c 



/7T 



(78) 



(79) 



Arbitrary additional terms corresponding to longitudinal fields may be added to the RHS's 
of (jT8"]) and ( 1791) without effect upon (164)) and (166)) ; here such terms are set to zero without 
loss of generality. 

The pair of equations (1761) and (ITT)) may be recast in the form of a single, generalized 
Hermitian eigenvalue problem for each value of Q in the domain (0, oo); viz., 



V x 



disp 



disp 



It 



V x 



j-\ee 



diss 



>,n) 
; r ,n) 

>,n) 



■pern . 

disp' 



r,fi) 



i - rr P (r, n) 



rs,(r,n) 





Don(r) 




n 


D Qn (r) 














B^n(r) 






B^n(r) 



-fr(r,n) iX(r,n) 
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D n „(r) 


j 




Bnn(r) 





(80) 



Mode field solutions and their corresponding, generally complex, eigenvalues Zq u are labeled 
by the index n. To obtain the normalization condition for the solutions to flHUl) . and thus 
the TP mode fields, we substitute the expansion (I64p into the ETCR (IB^|) . which we then 
evaluate using the expressions for the expansion coefficients (172"]) . (1751) . (178"]) and f[7§]) . with 
the result 

c Qn (t),cl, nl (t) 



427 r5(Q - SI') Z ^' Z ^ + ^ 

V ' IflTX 2 



X 



dK 



D 



f2n' 



B 



rS.(*.«) 

-fZ(r,0) C(r,n) 



In obtaining ([8"T]) we have exploited the identity 34] 
1 _ 1 _ 1 _ 1 



Q' - Q Q" - Q 



+ 



+ 



1 



B^n(r) 
1 



(81) 



Q - n» Q' - Q" 



(82) 



Comparison of (18 ID with the RHS of ( 1561) then implies 



n 



J Un\ 



X 



d^ 



D^(r) B^,(r) 



HL(r,fi) 



-pern . 
L diss 1 



■ diss 



>,0) C(r,n) 



D^n(r) 
Bcn(r) 



Ml <L 



(83) 



The procedure for constructing the D and B operators is then as follows. For each value of 
fi>0we solve (180]) with appropriate boundary conditions to obtain a discrete set of mode 
fields Do n (r) and B^ n (r), normalized according to (1831) . with eigenvalues Z^ n . Combined 
with the properties of the TP operators, these solutions then provide a complete description 
of the transverse electromagnetic field operators D and B through ( IBTT) and ( |62l . This 
description relies only upon the macroscopic susceptibility tensors and is therefore free of 
the ambiguities associated with the coupling tensors A^(r,f2). 



C. Other field operators 



In the modal polariton picture the TP mode fields may be employed to obtain explicit 
expressions for the noise operators. Specifically, by substituting (I72|) and (1731 into (163]) . we 
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may use the result to rewrite (13T|) and (|38|) . Comparison with (H3|) and (1441) in the Fourier 
domain then yields 



p( n )(r, U ) = " l) [rdiss(r^) • D wn (r) + • B wn (r)] c wn 

+^5Z5Z A A( r ' w ) -PAn(r,w) S Awn , (84) 

A n 

M<*>(r, W ) = - l) ra:(p,a;) ■ D_(r) + r™(r )W ) ■ B_(r)] c_ 

+^5Z£ A A( r ^)-PAn(r,w) s w (85) 



A n 



As alluded to in Section [III Bl these expressions clearly demonstrate how the deviation of the 
eigenvalue Zn n from in determines the contribution of each TP mode to the noise operators. 
We note that those parts of the noise polarization and magnetization operators associated 
with the TP operators only depend upon the macroscopic susceptibility tensors. 

The remaining electromagnetic field operators E and H may now be expressed in the form 
of mode expansions as follows. Substituting (|6ip . f )62|) . and the inverse Fourier transforms 
of flEU and (JS5J into and we may use (1551). (j75jl ; and (J7HD to write 



A n 



dfi 



dfi 



[A c A (r,fi) ■ p An (r,fi) s X nn(t) 

+A c A *(r,fi)-p An (r,fi) s\ nn (t) 



= E( ^ M) H^(r) + cL(t) H£(r) 



+Ar(r,n).pi n (r 1 n) 4 n „(t) 



(86) 



(87) 



Thus, in contrast to the transverse field operators D and B which may be completely 
described by the TP operators and their associated mode fields, the field operators E and 
H include contributions from the LP operators. 
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IV. CONCLUDING REMARKS 



We have presented a canonical quantization of macroscopic electrodynamics in a linear 
magneto-electric medium. The theory supports a wide class of magneto-electric responses 
characterized by the macroscopic susceptibility tensors T au (r,t), for a, v — e,m, which are 
Kramers-Kronig and Onsager consistent. The resultant electromagnetic field operators are 
expressed in a mode expansion representation and form a natural basis for the study of 
quantum optics in waveguide and cavity geometries involving dispersive and lossy magneto- 
electric media, while also paving the way for the inclusion of quantum optical nonlinearities 
in such structures. 

In practice, the available measured or modeled quantities are the susceptibility tensors 
T av '. Thus (H91) may be considered as a definition of the coupling tensors A^. However, the 
set of r"""s encompassed by our model is restricted. Specifically, the form of (|49|) implies an 
interdependency between the various susceptibility tensors due to the fact that the magneto- 
electric response is constructed from the same coupling tensors as the dielectric and magnetic 

n 

terms. Using the same methods as Horsley [27[ it may be shown that for any number of 
subsystems N, the following equality is satisfied: 

III- (r,c)| 2 < |r&(r, W )||r^(r,o;)|. (88) 

Thus, all media to which the present quantum theory applies must exhibit responses that 
satisfy the condition (ISSI) . 

As an example of the additional freedoms obtained by increasing the number of subsys- 
tems with which the medium is modeled, we first consider the case where N — 1. If the 
medium is isotropic we have 

rS.(r, w) = iJ {r cc (r, u)} 1 = ~A?(r, M) • Af (r, | w |), (89) 
IX(r, u ) = i J {r mm (r, u,)} 1 = ^A?(r, ■ A™(t, (90) 

where T ee and r mm are scalar functions and 1 is the 2nd rank tensor with Cartesian compo- 
nents 5ij. From ( 18"9"|) and (19T?|) acceptable forms of the coupling tensors are given to within 
an arbitrary unitary transformation by 

A^(r,fi) = [2^{P c (r,fi)}]M, (91) 
Af(r, n) = ±i \2JP {r mm (r, fi)}]' 1. (92) 
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The magneto-electric susceptibilities are then fixed by the dielectric and magnetic responses 
as 

IX s (r, co) = ±f-J? {P c (r, | W |)}' J {r mm (r, | w |)}* 1, (93) 
r^(r,c,) = -IX(r,u,). (94) 

This condition is relaxed by extending the model to N = 2 for which we have the general 
tensor relations 

I^Lfru,) = ~ [A;(r, ■ A^r, | W |) + A e 2 (r, | w |) ■ A c /(r, | W |)] , 
C(r, W ) = ~ [Af(r, |w|) ■ Ar(r, |w|) + A?(r, |w|) ■ A?* (r, |w|)] , 

r 4 ^ 

rSL(r, = ~ [Afr, | W |) ■ Ar(r, |w|) + A c 2 (r, | w |) ■ Af (r, |w|)] . (95) 

However, since all four coupling tensors contribute to all the susceptibility tensors in fl95|) . 
properties such as resonances which are present in the electric or magnetic susceptibilities 
must also manifest in the magneto-electric susceptibility. Now consider N = 3 in the special 
case where we set A 2 = A™ = to maintain the same number of coupling tensors in the 
model as before. Such a prescription yields 

rs.fr") = ~n ^ M) • A i>' M) + A 3( r > M) • A 3>> M)] > 
C(^) = t A "( r ' m) ■ A "*( r ' M) + ^fr M) ■ A 3>> M)] , 

rSLfr = ^A c 3 (r, | W |) • Ar(r, M). (96) 

Though the number of coupling tensors are the same as in the N = 2 case leading to ( |95l) . 
in fl96|) there are elements of the electric and magnetic susceptibilities that are constructed 
from A^ and A™ and thus decoupled from the magneto-electric susceptibility, which only 
involves A3 and A™. This allows for the inclusion of purely electric and magnetic effects 
which do not manifest in the magneto-electric response. 

Finally we recall that in all cases the coupling tensors are defined to within the unitary 
transformations Ua introduced in Section III BL Inspection of (1491) implies that the ambigu- 
ities associated with such transformations, however, have no effect upon the macroscopic 
susceptibilities which themselves are independent of XJ\. The noise operators defined in 
(jl5|) and (JIB"]) , and re-expressed in fl84l and (jSlj) . are likewise unaffected. More generally, we 
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expect that this independence with respect to Ua must apply to any physical result of the 
theory. 

In our quantization procedure we have employed the electric induction D and the mag- 
netic induction B as the canonical variables corresponding to the electromagnetic field. This 
prescription is equivalent to approaches involving the vector potential A and its conjugate 
momentum n, with D = IlandB = VxA (cf. the ETCR flU}). In relating the E and H 
field operators to D and B, the classical definitions of the polarization and magnetization 
are carried over into the quantum domain resulting in fl36l) . and the transversality of D and 
B is preserved. In this respect we are consistent with the work of Suttorp 13|. In contrast, 



Philbin 



ll| only retains the driven part of the polarization and magnetization in fl36|) . with 



the noise operators separated into additional source terms in the divergence equations for 
D and H; the operator D is therefore no longer transverse. 
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